The Absence of Cold Dust and the Mineralogy and Origin of the Warm Dust 

Encircling BD +20 307 

A. J. Weinberger 

Department of Terrestrial Magnetism, Carnegie Institution of Washington 
5241 Broad Branch Road NW, Washington, DC 20015 
weinberger@dtm. ciw.edu 
E. E. Becklin 

Department of Physics and Astronomy, University of California Los Angeles 
Box 951547, Los Angeles, CA 90095-1562 
becklin@astro.ucla.edu 
I. Song 

Department of Physics and Astronomy, The University of Georgia 
Athens, GA 30602-2451 
songOuga . edu 
B. Zuckerman 

Department of Physics and Astronomy, University of California Los Angeles 
Box 951547, Los Angeles, CA 90095-1562 
benOastro . ucla . edu 

ABSTRACT 

Spitzer Space Telescope photometry and spectroscopy of BD +20 307 show that all 
of the dust around this remarkable Gyr-old spectroscopic binary arises within 1 AU. 
No additional cold dust is needed to fit the infrared excess. Peaks in the 10 and 20 fim 
spectrum are well fit with small silicates that should be removed on a timescale of years 
from the system. This is the dustiest star known for its age, which is >1 Gyr. The 
dust cannot arise from a steady-state collisional cascade. A catastrophic collision of 
two rocky, planetary-scale bodies in the terrestrial zone is the most likely source for this 
warm dust because it does not require a reservoir of planetesimals in the outer system. 



1. Introduction 



A number of extreme collisions may have shaped the final states of the terrestrial planets. 
Impact by a Mars-sized progenitor on the early Earth, perhaps 30 Myr after the Solar System 
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formed, could hav e created the Mo on (e.g. ICanupl 12004 ). Other collisions may have stripped 
Mercu ry's mantle (IBenz et al.l [1988) or caused the hemispheric crustal thickness asymmetry of 
Mars (jMarinova et al.ll2008l ). Subsequent bombardment of the planets may not have generated 
planetary scale changes, but the "Late Heavy Bombardment," an era approximately 600 Myr after 
Solar System formation, either culmi nated a long sequence of bombardmen ts or resulted in a burst 
of impacts on the terrestrial planets (jGomes et alJl200a ; IStrom et al.ll2005l ). 



Evidence of cataclysmic impacts in other planetary systems is rare. Enhanced collisions should 
result in the production of copious small particles, i.e. dust, which can be observed by their re- 



radiation of stellar light (jKenvon fc Bromlevl 12004 ) . Stars older than 100 Myr with substantial 



Mever et al. 


2008: 


Trilling et al. 


2008: 



2009 



Koerner et alJ 120101 ). 



At any given age, stars that are dustier than average could be interpreted as having had a recent, 
stochastic, event, but many could also be the result of differing initial masses and configurations of 
planetesimals. Only the extremely dusty disks at any given age must be from giant collisions. 



BD+20 307 is a relatively anonymous field star at 96^q pc ( van Leeuwenl 120071 ) from the Sun. 



It is actually a tidally locked spectroscopic binary wit h a period of 3.5 d composed of two n early 
identical late F-t y pe st ars with a mass ratio of 0.91 ([Weinberger! 120081 ). IWeinbergerl (J2008J) and 
Zuckerman et alJ (120081 ) investigated the stellar properties, in particular the weak and disparate 
Lithium abundances in the two st ars, and estimate d an age > 1 Gyr. At that age, its dust opacity, 
estimated by Ltr/L* to be 0.033 (ISong et al.ll2005l ). is orders of magnitude larger than every other 
star of its age. In fact, such dustiness would be extraordinary even for a star 10-100 times younger. 

A sma l l clas s of very (wa r m) d usty objects now exists in the literature (see summaries in 
Rhee et alJ (|2008l ) : iMoor et alJ (|2009)). The extreme example s of these and also some disks wit h 
strange infrared spectra such as HD 172555 ( Lisse et al. 2009 ) and 51 Oph; ( Ancker et al.l l200ll ) . 
have been attributed to impacts of planet-sized bodies. 

Here, we investigate the unusually dusty system around the star BD+20 307 using the Spitzer 
Space Telescope. These observations constrain the size distribution and location of a very large 
amount of dust and provide perhaps the best evidence yet for large-scale collisions in another 
system. 



2. Observations and Data Reduction 



We observed BD +20 307 with Spitzer with all three instruments - IRAC, MIPS, and IRS. 
An observing log is given in Table HJ 



2.1. IRAC 

IRAC data were taken in subarray mode, which yields 64 individual images. The S14 pipeline 
reduction was used, and the "Basic Calibrated Data" (BCD) images were mosaiced with the Spitzer 
Science Center's MOPEX software. There are two photometric corrections to consider. The first 
is the "Pixel phase dependent photometric correction" that corrects the pixel response function for 
where the peak response is within a pixel. This correction was applied to channel 1, and amounted 
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to 1.7%. The second is the array location dependent photometric correction, due to the fact that 
IRAC is flat-fielded on the Zodiacal spectrum. From the overall shape of the BD +20 307 spectral 
energy distribution (SED), we know that its spectrum is shallower in slope than Rayleigh- Jeans 
from 3 -7 um but not as flat as the Zodiacal emission. The correction was applied at Channel 1, 
and amounted to 0.7%. It was also applied to channels 2 and 3, but made less of an effect, and 
was not applied for channel 4. Aperture photometry was performed in the radius 2, 3, and 5 pixel 
radii apertures given in the IRAC handbook, and the standard deviation of these taken as the 
measurement uncertainty. 

Finally, the IRAC data need to be color corrected from the default vF v = constant calibration 
to the true BD+20 307 spectrum. This effect is small for Bands 1 and 2, where the source spectrum 
is approximately a 1000 K blackbody, and the correction is 0.7%. For bands 3 and 4, we integrated 
the observed IRS spectrum over the IRAC filter bandpasses. This yields a correction of 3% at Band 
3 and 71% for Band 4. The Band 4 correction is unreasonably large; the difference between the 
integral under the IRS spectrum and the IRAC Band 4 (un-color corrected) flux density is only 
26%. When integrated, the very large spectral slope at 8 //m due to the silicate feature that peaks 
at ~10 jim, causes a large change in color correction for any small change in spectral transmission 
response of the system (telescope + filter, etc.). This transmission does not appear to be well 
enough known to provide an accurate color correction; therefore, the Band 4 flux was not used for 
further analysis. 



2.2. MIPS 



At 24 um, the pipeline produced post-BCD mosai c was used for aperture photometry using a 6 
pixel aperture and corrections given in lSu et al.l (120061 ). At 70 am, MOPEX was used to mosaic the 
pipel ine filtered BCDs. Ap erture photometry was performed with the 35" aperture and corrections 
from [Gordon et al.l (|2007l ). At 160 um, the pipeline mosaic was used. Because at this band no 



Table 1. Summary of Observations and Photometry 



Instrument 


Band 


Date 


Integr. time (s) 


Photometry (mJy) 


IRAC 


1 (3.5//m) 


2005-08-20 


0.64 


285.0 ± 7.3 


IRAC 


2 (4.5/zm) 


2005-08-20 


0.64 


210.4 ± 4.8 


IRAC 


3 (5.7/xm) 


2005-08-20 


0.64 


165.3 ± 12.0 


IRAC 


4 (7.9/«n) 


2005-08-20 


0.64 


410.2 ± 22. 5 a 


MIPS 


1 (24/zm) 


2007-01-21 


2.62 


441.0 ± 0.8 


MIPS 


2 (70/xm) 


2007-01-21 


167.8 


28.6 ± 1.9 


MIPS 


3 (160/zm) 


2007-01-21 


209.8 


<22 


IRS 


(5.2-14.5/im) 


2006-01-15 


12.6 




IRS 


1 (9.9-19.6/im) 


2006-01-15 


188.7 




IRS 


3 (18.7-37.2/xm) 


2006-01-15 


1219 





Not color corrected 



- 4 - 



source is seen at the nominal location or elsewhere in the field, an upper limit to the source flux 
was placed by doing a raster of aperture positions on the image around the nominal location and 
calculating three times the standard deviation of these. 



2.3. IRS 

For all channels, the pipeline version 14 data were used. The short-low (5.2 - 14.5 /jm) BCD 
data were cleaned with IRS_CLEAN software to remove bad pixels. Each cleaned BCD was sky 
subtracted with uncertainties propagated in quadrature. The spectra were extracted with the 
Spitzer Science Center's SPICE software using the standard extraction parameters. 

The short-high (9.9 - 19.6 /tm) BCD data were cleaned with IRS_CLEAN and the unaltered 
campaign pixel mask. The resulting cleaned BCDs were medianed to make a final spectrum that 
was extracted with SPICE and the standard extraction parameters. The two nods were averaged 
together. Orders were trimmed to eliminate as many points with undefined uncertainties as possible. 

For the long-high (18.7 - 37.2 /mi) data, separate observations on nearby sky were taken to 
help in elimination of rogue pixels. A median sky frame was created with the standard deviation of 
the 20 sky files as its uncertainty. This frame was then subtracted from each of the object spectra, 
with uncertainties propagated. 

The campaign rogue mask was augmented by any pixels that were flagged in 15/20 of the 
individual spectra. The combined mask was then used to clean the sky subtracted object spectra. 
The final cleaned BCDs for each nod were median combined and then extracted with SPICE. The 
spectra from the two nods, which agree to within their statistical uncertainties, were averaged. 

The spectra from all channels were combined by using their overlap regions - short-high to 
short-low spectrum in the region 11.5-14 /jm and long-high to short-high at 19.02-19.47 /mi. The 
complete spectrum was then normalized to the MIPS channel 1 photometry by integrating the 
spectrum under its transmission curve. A check that this process works well is that integrating 
the final normalized spectrum under the IRAC channel 3 filter yields 155 mJy compared with the 
measured 165 ± 12 mJy from the photometry reported in Table [TJ 



Results 



Results of the new Spitzer photometry are given in Table [TJ As first reported in ISong et al 



(|2005l ). the spectral energy distribution of BD+20 307 at wavelengths longer than 5 /mi is dominated 
by emission by small silicate particles. We show the overall SED including the new Spitzer spectrum 
and photometry along with previous ground-based photometry in Figure HJ 



3.1. Photosphere Modeling 

We model the combined flux density of the nearly identical central stars with a single Kurucz 
model atmosphere fit to the Tycho-2 and 2MASS photometry at 0.4 - 2.2 /mi. The best fit model 
has a gravity of log(g)=4.5, T e g=6000, and, at the parallactic distance of 96 pc, a luminosity of 
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1.94 L/0. IZuckerman et al.l (j2008l ) finds the luminosity ratio of the two stars as 0.78, which would 



imply individual luminosities of 1.10 and 0.84 L . Table [2] gives the amount of photosphere in each 
Spitzer band and the remaining disk flux. The model photosphere was also subtracted from the 
IRS spectrum. 



3.2. Spectral Fitting 

We proceed by assuming a simple model for the composition and location of the circumstel- 
lar dust and asking whether it is a plausible match to the data. The presence of both broad 
and sharp spectral feat ures indicates emission by both amorphous and crystalline materials (e.g. 



Campins fc Rvanlll989l ). 



The shortest observed peak, at 9.7 /im, is matched well by the wavelength of the peaks in 
measurements of glassy silicates, although its width is narrower than that produced by even very 
small (< 0.1 /mi) grains. The breadth of the 16-20 /mi peak is better fit by larger, 1 /mi grains. 
We chose to set the minimum grain size to 0.5 /mi, which is approximately the blow-out size for 
radiation pressure from the stars. We then use a power law distribution of sizes rising from 0.5 /mi 
with slope -3.5, but the fit does not materially change whether the power law distribution or just 
the minimum-sized grains are used. Smaller grains result in narrower peaks that improve the fit at 
10 /mi but make it worse at 20 /mi. 

We take the simple model to be dust at a single distance from the star and composed of four 
types of particulates: small amorphous dust of olivine composition (MgFeSiO^ hereafter called 
amorphous olivine), small crystalline olivine (forsterite; Mg2Si04)), small amorphous pyroxenes 
(Mgo.5Feo.5Si03), and large grains (i.e. those that behave as blackbodies) . These are common 



constituents used to model solar system comets and circumstellar disks (e.g lLi et al.ll2004l ). 



All small particles (size << wavelength of observation) have absorption efficiencies that fall 
faster than Rayleigh-Jeans (i.e. oc A -2 ) at wavelengths much larger than the particle size. Therefore, 
some bigger grains are necessary to match the measured flux densities at 20-70 /mi. Again, as a 
matter of taste, we have chosen to keep the grain size distributions for the olivines and pyroxenes 
the same, as that seems the most likely physical outcome of a collisional cascade, and introduce a 
blackbody component to the fit. The minimum size of the pyroxenes is not well-constrained by the 



Table 2. Photosphere subtracted flux densities and model fits in mJy 



Band Photosphere Disk flux Model (mJy) 



IRAC1 


261.0 


24.0 


9.3 


IRAC2 


167.0 


43.4 


28.8 


IRAC3 


107.1 


58.2 


70.8 


IRAC4 


60.7 


349.5 


293.2 


MIPS1 


6.8 


434.2 


424.4 


MIPS 2 


0.8 


27.8 


40.2 


MIPS3 


0.3 


<22 


8.7 
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shape of the spectrum, because pyroxenes have relatively broad features. A larger minimum size 
for the amorphous pyroxenes could be traded off against the blackbody grain component. 

We assume that the disk dust is optically thin, so that each dust grain is in thermal equilibrium 
based on its own absorption efficiency and that the various emission components simply add. The 
Spitzer 3.5, 4.5, 5.7, and 70 pm photometry and the combined 5.6 - 25 ^m spectrum are fit 
simultaneously. 

We write the emission due to the disk dust in the optically thin case as: 



F u = T,AiKiB Uti (d) 



(1) 



where the sum is over i independent dust species each with their own mass absorption coeffi- 
cients, Ki = 3Q a bs/4:ap where Q a bs is the absorption efficiency, a is the grain radius, and p is the 
grain density. The dust grains are not assumed to be in thermal contact; thus each species has 
its own temperature determined by its Q a bsi distance from the star, and the luminosity of the star 
by assuming it is in thermal equilibrium. These temperatures are calculated self-consistently by 
making distance from the star the free parameter. The other free parameters are the number or 
surface density of dust grains of each type, i.e. the A4. 

A major choice is what absorption efficiencies to use. Various authors have taken different 
approaches to this problem. Because we do not leave temperature as the free parameter, we 
need genuine ultraviolet and visible absorption efficiencies, as this is where the star emits and the 
grains absorb most efficiently, in order to calculate temperatures. For this reason, we opt for Mie 
calculations. 



Mass absorption effic iencies from the group at Kyoto Pharmaceutical University (jKoike et al 
2000; ISogawa et al.l 120061 . e.g.) are available only in the infrared. Optical constants from the 



Jena group are provide d into the UV- visible, but are only accurate in the near- infrared and longer 
(Dorsch ner et al.l Il995l . e.g.). For the amorphous olivine, we adopted the optical constants of 



Li et al 



(120041 ) which combine the Jena measurements with ultraviolet data. We then used a Mie 



code 1 to calculate the mass absorption opacities. 

For crystalline olivine, shape is a major determinant of spectrum. The Jena group has measured 
the triaxial indices of refraction, but how one uses these to calculate the absorption i s a matter of 



taste . Assumption of spherical particles results in sharp peaks that are not observed (IFabian et al 



2001 



The Jena group provides two computations of continuous distributions of ellipsoids (CDE) 
in the Rayleigh limit. We combined these with the ultraviolet efficiencies calculated in the Mie 



limit from the constants of lLi et al.1 (J2004J). We tested both the CDE1 (weighted toward spherical 
grains) and CDE2 (all shapes weighted equally) calculations against our spectra and found that 
CDE1 was a bett er match. It is also CDE1 that fairly closely matches a real grain measured in the 



Jena laboratory ([Fabian et al.l 12001 



It is because of the flexibility in spectral shape granted by choosing the exact distribution of 
grain shapes, that we have not attempted a multi-component fit that would reduce our chi-square 
to a level approximately equal to one. We believe that the large number of free parameters to 



http://www-atm.physics.ox.ac.uk/code/mie/index.html 
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choose from makes this exercise of limited value. 

To find the best fit, we compare the model excess as described above to the photometric (2.16 
- 70 /jm) points and IRS spectrum in the classic \ 2 formalism. We first perform a search over a 
large ran ge of parameter s pace for the free parameters using the Amoeba routine from Numerical 



Recipes (IPress et al.lll992l ) as implemented in IDL. Once the best fit has been found, we then loop 
over a small range of values for the free parameters to determine the shape of the x 2 an d determine 
a better minimum. 

The best fit is shown in Figure El The best distance to the dust is 0.85 AU from the star. 
The calculated temperatures of the grains are 461 K for the pyroxene, and 460 and 507 K for the 
amorphous and crystalline olivine silicates, respectively. The blackbody temperature is 358 K. The 
reduced chi-square value of this fit is 19.7. Obviously, the reduced chi-square is not near one, but 
nonetheless, the fit is capable of explaining all of the main features of the spectrum including the 
excess emission at 3-5 ^m the whole of the 70 //m emission, and the locations of peaks ((10, 12, 17, 
19, and 24 fim) in the 7-35 fim spectrum. The overall match between spectrum and the model is 
good except in the heights of the peaks at 10 and 17 //m and the continuum at 25-30 fim. For the 
issues of most relevance, namely the temperature of the dust, the presence of small grains, and the 
crystallinity fraction, this model adequately represents the spectrum. Its predictions in the various 
filters are given in Table [21 The best fit parameters are given in Table El 

The contributions by mass are 50% amorphous olivines, 30% amorphous pyroxenes, and 20% 
crystalline olivines. The actual masses in these small grains range from 0.7 to 1.7 x 10 22 g. 
Obviously, most of the mass would be in the blackbody grains. A minimum mass for these can be 
set by setting the minimum grain radius for one of the above grains to look like a blackbody, or 
~30 /im. Then, the minimum mass is 6 x 10 23 g. 



4. Discussion 

The remarkable results of the Spitzer data are the lack of cold dust - all of the 70 ^m emission 
arises from the same close-in population of dust producing the bright 8-30 /im silicate features, 
and that the dust that does orbit the star is extremely close-in - within 1 AU. It is the most dusty 



Table 3. Best Fit Model 



colheadComponent 


A, 


T(K) 


Ltr/L* 


Amor. Olivine 


1.0 a 


460 


0.010 


Amor. Pyroxene 


0.653 a 


462 


0.006 


Crys. Olivine 


0.821 a 


508 


0.004 


Blackbody 


1.64e-4 


358 


0.011 


distance 


0.85 AU 







a These coefficients, defined in Section [3T2| are nor- 
malized such that Amor, olivine is 1.0. 
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star for its age known and the one with the warmest dust. 



There are now about ten known "warm" debris disks around FGK stars ( Rhee et al 



Moor et al 



2009 



Fuiiwara et al 



2010 



2008 



Melis et al.ll2010i ). Most of these have dust temperatures of 



~150-200 K and Lj#/L*<5x 10~ 4 . For stars older than 1 Gyr, BD +20 307 has more than two 
orders of magnitude more dust opacity than the others. The closest star in terms of dustiness is 
the much younger Pleiades member HD 23514. 

The timescale over which the BD +20 307 disk shou ld evolve are short. The treatment of 
radiation pressure on small silicates in iBurns et al.l (|1979l ) and scaled to the slightly larger mass 
and luminosity of BD+20 307 suggests that all grains smaller than 0.5 //m will be blown away on 
an orbital timescale, which at 0.85 AU, assuming a central mass of 2.2 M©, is 0.5 yr. 

Taking the model fit L/^/L* = 0.032 as the surface density of the grains yields a collision time 
of 2.4 yr. This is only a few times longer than the orbital t ime, so the small gr ains observed in the 
spectrum will be created and then lost nearly immediately. IWvatt et al.l (|2007l ) argue that massive 
planetsimal belts grind themselves down by collisions quickly, so it is not possible to maintain very 
large amounts of dust in the terrestial planet region of old stars. They included BD+20 307 in the 
list of stars for which it was impossible to generate the present dust through collisional grinding, 
even using the IRAS 60 /um upper limit to constrain the cold planetesimal population and for 
an age of 300 Myr. At the current best estimate of its age of at least 1000 Myr , its m aximum 
Ltr /L* in steady state (see their Eq. 20) is 7.4 x 10 -8 . Recently. iHeng fc Tremaind (2010j) confirm 
the lWvatt et al.l (|2007l ) result that BD +20 307 cannot be from a steady-state collisional cascade, 
although their models allow most of the other warm disk sources to come from a warm planetesimal 
cascade. 

We do not detect any cold dust; our best-fit model for the hot dust actually overpredicts the 
measured 70 //m flux density and is slightly higher than the 1 a upper limit at 160 /im (Figure 
[3j. A check on the reasonableness of the lack of cold dust, is that if we take the average flux 
density in the IRS spectrum of 0.4 Jy and extrapolate this as a Rayleigh-Jeans fall-off from the 
mean wavelength of the IRS spectrum, 20 /im, we would predict a 70 mum flux density of 0.032 
Jy, which is again higher than measured. 

We can place a conservative upper limit on the amount of cold dust by finding the luminosity 
of a dust ring whose spectrum peaks at 160 /im (i.e. temperature of 33 K) and produces no more 
than a 10% increase in the flux density at 70 /im - such a cold disk would produce an additional 
6 mJy at 160 /im, which would still fall under the measured 3 a upper limit and would have a 
L/r/L* = 3 x 10 -5 . Therefore, in the context of the Wyatt et al. model for planetesimal belt 
self-grinding, the lack of detectable cold dust rules out any belt within 200 AU that is losing mass 
at a rate sufficient to account for the hot dust. 

Late heavy bo mbardment models per se are not attractive analogies for BD+20 307. In the 
(jGomes et al.ll2005l ) scenario, a late crossing of the 2:1 resonance between Jupiter and Saturn swept 
resonances across the asteroid and Kuiper belts and dislodged many small bodies that then impacted 
the terrestrial planets. For a similar scenario to work for BD+20 307, the resonance crossing must 
be delayed even longer, until 1 Gyr, but that is likely achievable by tuning the initial configuration 
of planets. More problematic is that the Kuiper Belt prior to the late heavy bombardment, which 
was massive enoug h to generate giant planet migration, would itself produce significant cold dust 
(|Mever et al.lboOTl ). 
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The remaining hypothesis for the source of the hot dust around BD+20 307 is a giant impact. 
Detailed models exist of one such impact - the Moon- forming event. About a quarter of the 
material that forms the circum-terrestrial disk out of which the Moo n formed sta rts very hot, at 
temperatures >1000 K, or high enough to vaporize silicates and metals ()Canupll2004l ). Furthermore, 
a few percent or more of a Lunar mass can escape the Earth-Moon system. The small particles 
around BD+20 307 are ~20% crystalline. This level of crystallinity is well within range of T Tauri 
star dust; in contrast, differentiated bodies such as large asteroids or planets should be nearly 100% 
crystalline. That the observed dust is largely amorphous suggests that the impact heated the bulk 
of the material to temperatures ~1000 K and that any recondensation took place under low vapor 
pressure conditions where the molten silicate cooled quickly and therefore did not re-anneal into 
crystalline form. 

A te rrestrial planet sy stem may be stable over a Gyr but not permanent. The Planet V 
models of IChambersI (|20071 ) show that in a configuration such as the Solar System's, a fifth planet 
can survive in the inner Solar System for a Gyr and create collisions of terrestrial planets at many 
hundreds of millions of years after the system formed. The configuration of the exoplanet system 
will determine how long it lasts; small planets further from the star take longer to go unstable 
(Chambers, personal communication). In the case where a tenuous asteroid belt still remains in a 
planetary system, the trajectory of a rogue planet can destabilize the belt and produce collisions 
between planetesimals or bombardment of the planets. Over several Gyr, o ur current terrestrial 
plane t system may be unstable and result in collisions between the planets (jLaskar fe Gastineau 
200fll ). 



Recently, iMatranga et al.l (|2010i ) suggested that the changing separation of close binaries due 
to angular momentum loss could sweep resonances across planetary systems and dest abilize their 
planet esimal belts. The stars in BD+20 307 are separated by 0.06 AU at present IWeinberger 
(|2008l ). and thus the disk radius is 14 x the binary separation. They have likely been in a tight 
circular orbit for billions of years IWeinbergerl (|2008l ) and their relatively low X-ray activity likely 
makes their angular momentum loss relatively slow. 

How common are giant collisions amongst sun-like stars? We can estimate the number of giant 
impacts per star N g = f^AL -1 where /* is the fraction of stars observed to have hot dust, L is the 
lifetime of the collision products, and A is the age of the stars surveyed. 

To estimate the lifetime of the dust, we invert the rtWyatt et al ] l2007h expression for the amount 
of dust that can be produced by a collisional cascade over time (their equation 20) to solve for the 
lifetime of a collision that starts in a giant impact. The most significant unknown is the internal 
strength of the large planetesimal(s) involved, but the binding energy is likely closer to that of a 
differentiated planet than a pure rubble pile, i.e. more like 2xl0 5 than 200 J kg" 1 . Therefore, 
a planet-sized impact could generate a cascade that lasts up to 80,000 yr. This result is similar 
to the 100,000 yr lifetime for a collisional cascade started by planetary embryos as calculated by 
(jMelis et alJboiol ). 



Surveys with Spitzer have looked at all Sun- like stars within 25 pc lKoerner et al.l (|201Cj), for a 
sample of about 800. T he average age of th is population is >3 Gyr but it includes a small fraction 
of much younger stars (jWright et al.ll2004l ). While there are five stars with warm dust in this 
sample, we note that these have 24 /mi flux densities <10 times their photospheres and Lir/L* < 
10~ 3 , and so are not nearly as extreme as BD+20 307. Moreover, their average age is <1 Gyr. 
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Their warm dust may either be from young planetesimal belts (e.g. j3 Pic, which is in this group) 
or from collisions fed by outer belts and constrained by planets (e.g. r/ Corvi; Wyatt et al. 2007). 

IRAS was capable of detecting every A-G star with as much hot dust as BD+20 307 out to 
~280 pc, although it is possible some of these would not have been noticed. There are ~600,000 
stars stars of earlier spectral type than K0 in that volume. Four stars (including BD +20 307) 
have F24/F* > 10. This fraction of host dust stars implies a giant impact rate greater than 
0.2 impacts/star during its main sequence lifetime, i.e. implies that such impacts are common. 
Obviously, these are small number statistics. The WISE mission currently surveying the sky at 
mid-infrared wave lengths will gre atly improve the statistics. Its sensitivity is 500 times better than 
IRAS at ~24 //m (|Liu et al.ll2008l ). and so it should provide a truly complete census of isolated hot 
dust sources out to a few kpc. 

It is difficult to search for planets around BD+20 307 with the radial velocity technique because 
of the spectroscopic binary. However, if planets do encircle this star, they may trap small bodies in 
resonance. An enhancement in the planetesimal density might increase the liklihood of a destructive 
collision such as the one observed. The star HD 69830, also an IRAS source, and > 1 Gyr old, was 
shown with Spitzer to harbor a ring ver y similar in size and t emperature to that around BD+20 
307, but containing 150 times less dust (jBeichman et al.ll2005l ). A system of three Neptune mass 
planets was then found to encircle HD 69830, wi th the dust eithe r residing between the 2nd and 
3rd planets from the star or exterior to all three (jLovis et al.ll2006l ). 



Whatever the scenario that orginated the collision in the BD+20 307 system and others, it 
is interesting to note that the spectra of other warm dust systems are no t identical. HD 69830 , 
for example, has much more prominent 20 fim crsytalline silicate peaks ( Beichman et al. 2005) 



and HD 23514 has an unusal broad spectral feature at 9 rather than 10 //m (jRhee et al.l 120081 ) . 
HD 172555's narrow 8 /im peak and featureless 20/xm regions have been attributed to SiO. These 
various compositions may result from a combination of progenitor composition, which would affect 
the primordial abundances of different types of silicates, and geometry of the impact, which can 
change the melt fraction and condensation properties of its debris. 



5. Conclusions 



We have detected a large amount of small silicate dust at 0.85 AU from the >Gyr old star 
BD+20 307 but no cold dust that could be generated tens or hundreds of AU from the star. A 
catastrophic collision of two rocky, planetary-scale bodies in the terrestrial zone is the most likely 
source for this warm dust because it does not require a reservoir of planetesimals in the outer 
system. Furthermore, the high amorphous silicate content of the dust suggests that the impact 
must have caused a high melt fraction and amorphization of what was likely a largely crystalline 
parent body. 



This work is based on observations made with the Spitzer Space Telescope, which is operated 
by the Jet Propulsion Laboratory, California Institute of Technology under a contract with NASA. 
Support for this work was provided by NASA through an award issued by JPL/Caltech and a 
contract to Ames Research Center for the SOFIA program. We acknowledge support from NASA's 
Astrobiology Institute to the UCLA and CIW nodes. John Chambers and Sarah Stewart provided 
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Fig. 2. — Best fit model to the Spitzer spectrum (black line), derived as described in the text. 
The various dust components (purple: amorphous olivine, cyan: crystalline olivine, dark blue: 
amorphous pyroxene, and red: blackbody grains) sum to make the total spectrum (gold line). All 
components are situated at the best fit distance of 0.85 AU. 
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Fig. 3. — Spectral energy distribution of BD+20 307 showing the visual and near-infrared photom- 
etry from the Tycho-2 and 2MASS catalogs along with the Spitzer photometry and spectroscopy 
from Figure 1 and our best-fit model from Figure 2. We do not detect any cold dust; our best-fit 
model for the hot dust actually overpredicts the measured 70 /im flux density and is slightly higher 
than the 1 a upper limit at 160 //m (the 3 a upper limit is plotted). The model does a reasonable 
job of reproducing the data for the issues we care about, namely, the dust location (temperature), 
the presence of small grains, and the crystallinity fraction. 
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